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Motivation:  Non- 
Station  ary  Interference 


•  Rapidly  changing  clutter  locus  with  a  circular  array 
or  bistatic  radar  system 

•  Presence  of  hot  clutter  due  to  an  airborne  jammer 

•  Use  model  of  non-stationary  interference  to  derive 
new  filter 

•  Use  small  sample  support  to  reduce  effect  of 
non-stationary  interference 
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Data  Model 


•  M  antennas,  N  pulses 

•  Target  in  primary  range  bin  p 


xp(t)=  ba(e  )eM  +  cp(t),  t=  0,1,2, 

_ / 


spatial  steering  vector 


clutter ,  jammer,  noise,  etc. 


Space-Time  Slice 

Xp  =  [xp(0)  xp(1)  xp(N 
=  i ba(e  )vT(<o )  +  Cp 


-1)] 


temporal 
steering  vector 


=  [l  e 


je»/Ts 


J  (N  -1)®/TS 
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Data  Model  (cont.) 

•  Vectorized  Forms 

1.  xp  =  vec(Xp ) 

=  bv(&)®  a(e)  +  cp 

2.  xp  =  vec (Xp ) 

=  ba(e)®v((o)  +  cp 

•  Secondary  Data 

{ck}  k  =  1,...  ,NS  k*p 
E (ck )  =  0  ,  E(ckc*k)  =  R 

_ / 


interference 

covariance 


target-free 
range  bins 
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Space-Time 

Autoregressive  Modeling 


Define 


L-l 


H(z1)=£H,.z 


i=  0 


Model:  for  some  L, 


M  ‘  x  M  matrices 


H(z  1  )ct (0  -  H0c* (0  +  HjC* (r - 1)  + ...  +  H,  , C;  (/ -  L+  1) 
=  £*  (0 

is  spatially  and  temporally  white 
•To  estimate  H(z_1),  solve 

NS  N 

Hmin  E  I  HuX 

H°’  ’Hl  *=1  i=L 
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STAR  filter  attempts  to  minimize  clutter  power: 


dimension 
M'(N-L+1)  x 


HL  l  hL  2 


q  |  Ck  -  M-Ck 


HL- 1  hl 


Span(J{)  orthogonal  to  clutter  subspace  if  it  dominates  white  noise: 

R  =  H  XQH  '*+  0  2I 

clutter  &  jamming^  '  white  (thermal)  noise 

so  we  project  onto  the  orthogonal  subspace  using  a 
matched  subspace  filter: 

x'  =  H  *(H  H  J1  H  x 


banded  block  Toeplitz 
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Algorithm  Summary 


1 .  Use  SVD  on  secondary  data  to  solve  for  [H0  H,  HL_ ,  ] 

computational  order:  0(NsM2L2 (N-  L+ 1)) 

2.  Form  Jf  and  filter  data:  P  ,x 

ji  p 

computational  order:  0{M'M2L2(N-L  + 1)) 


3.  Perform  regular  beam  and  Doppler  filtering  for  detection 

computational  order:  negligible 


Resultant  test  statistic  is  (v®  a)  P^.xp  not  (v®  a)  R-  xp 
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Prior  Work 


•  Vector  AR  models  used  previously  for  clutter 
modeling  by  Michels,  Rangaswamy,  etc. 


•  Standard  STAP  filters  extended  to  handle 
range-varying  and  hot  clutter  models  by 
Zatman,  Rabideau,  etc. 


•  Matched  subspace  detectors  used  for 
subspace  interference  by  Scharf 

•  Here,  we  extend  the  parametric  model  to 
^  handle  the  non-stationary  interference 
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To  improve  performance  at  short  ranges,  use  linearly 


•  Analogous  to  ESMI  technique  of  Hayward 

•  To  normalized  the  noise  subspace 

I 

V  (Ns  +  2)(A/S  +  1) 

_ _ _ / 
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Range-Varying  STAR  Filter 


•  Minimize  clutter  power  assuming  linearly  varying 
statistics 


ek  ~  ^ 


u k 

akc 


where 


J-C  from  STAR  filter 
T? 


Extended  STAR  filter 
coefficients  AJ( 


jS 


•  Filter  data  with  matched  subspace  filter 


k=0  for  primary 
range  bin 
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Define 


r  cjl+i) 


Ck(  1) 


Estimate  filter  coefficients: 


ck(N)  1 


ck(N-L)\ 


[H0  ■■■  Hl_ j  a H0  ■■■  aHl_j  ] 

as  the  left  singular  vectors  with  the  M’  smallest 
singular  values  of  the  extended  data  matrix: 

r  C-NJ  2  CNJ2  1 

CC  N  s  _  Ot  Ns  _ 

_  W  M  /  o  ••  W  Af  /  o 
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Normalized  Doppler 


I3TO  ECEn 


20  element  circular  array,  18  pulses 
SCR  =  -58dB  SNR  =  10dB 


Primary  data  vector  snapshot  at  20  km 


-SO  -SO  -40  -20  0  20  W  60  30 

Azimuth: 

_ 


4  tap  ESTAR  filter, 

20  secondary  snapshots 


40  40  -40  -20  c  20  40  $0  _ 60 

Azimuth 
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Computational  Comparison 


Some  typical  numbers:  M  =  20,  N  =  18,  M’  =  20 

•  STAR  Filter  (l_=5):  0(1 40,000A/S )  +  0(2,800,000) 

«  ESTAR  Filter  (L=41: 

0(4NSM2L2(N-  L+l))+  0(M  M2L2(N-L  +  l)) 
=  O(384,000A/s)+  0(1,920,000) 

•  Extended  PRI  staggered  algorithm: 
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Average  SINR  loss 
is  the  area  between 
the  curves 
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V- 

2 


20  km  range 


6  7  E 

Filter  order  {L> 


10  II  M 


30  km  range 


Ns=50  training  vectors  -  2  km  training  window 
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Average  SINR  toss  (dB) 
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Performance  with  Range-Varying  Weights 


20  km  range  30  km  range 


L=5  for  STAR  filter  -  L=4  for  ESTAR  filter 
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3-D  STAR  Filter  for 
Hot  Clutter 


•  Update  filter  for  each  new  pulse  received 

•  Derive  slow-time  varying  STAR  filter 

•  Can  be  used  with  intrinsic  clutter  motion 

•  Add  fast-time  matrix  taps  to  exploit  correlations 
across  range  bins 

•  Additional  filter  taps  help  mitigate  mainbeam 
jamming  signals 
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Slow  Time-Varying 
STAR  Filter 


•  Same  structure  as  the  STAR  filter  but  with  new 
coefficients  for  each  pulse 


•  Additional  sample  support  required  due  to  additional 
parameters  to  model  slow-time  variation 
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3D-STAR  Filter 


•  Use  slow-time  varying  STAR  filter  to  model  correlation 
across  pulses 

•  For  some  fast-time  filter  order  J, 
model  the  fast-time  correlation  as: 

j-i 


k  =  J  + 1,  ,P 

. _ 1 


number  of  fast-time  samples 


subscript  denotes  which  fast-time 
sample  J{  is  associated^  with 


used  to  whiten  data 


•  Similar  to  a  2-D  vector  AR  model  with  the  slow-time 
taps  changing  with  each  pulse 


J 
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Estimation  of  Parameters 


•  Define 


9k(t)=\ 

1 

+ 

o 

<5,(0  = 

rgfc+J-i(0 

1 

L  c*(f)  J 

9k+P-j(  0J 

Least  squares  solution: 


/v. 


min£  H(t)Gk(t) 

H(t '  fr=l 


•  New  minimization  for  each  slow-time  step 
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Filtering  the  Primary  Data 


•  3D-STAR  filter  can  be  written  as: 


•  Project  out  the  interference  using  3D  matched 
subspace  filter 


r  x‘  i 

X 

X 

p 

*  ,  * 

=  H  (HH  M 

p 

=  pr 

p 

*P-p+ 1 

_*p-p+ 1_ 

_XP-P+ 1_ 

Highly  structured  nature  of  subspace,  small  sample 
support  make  full  3-D  STAR  solution  feasible 
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Computational  Comparison 

Some  typical  numbers:  M  =  20,  N  =  18,  M’  =  20,  P  =  3 

•  STAR  Filter  (L=7):  0(235, 000A/S )+  0(4,700,000) 

•  3D-STAR  Filter  (L=2):  0{N S{MLJ)2 {N -  L  +  \\P -  J  +  \)) 

+  0{M'  (MLJf  0 N-L+ 1  )(P  -  J  + 1)) 

•  >2:  0(2 1 8,000 Ns)+  0(4,350,000) 

•  J5I:  O(82,000A/s )+  0(1,630,000) 


•  Optimized  3D-post-Doppler  algorithm: 

0{Ns(MKP)\N- K+  1))+  0{r\(MKPf(N- K+  1))=  0(5 1 8,000A/S )  +  0(70,000,000) 


#  of  sub-CPIs  =  3 


rank  of  sub-CPI  covariance  =  135 
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12- 


Direct  path  jamming  signal 
in  mainbeam:  JDOA=1° 


Only  multipath  component 
in  mainbeam:  JDOA=-20° 


Ns=100  training  vectors  -  4  km  training  window 
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Average  SINR  toss  (dB) 
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Hot  Clutter  Examples 

Direct  path  JDOA=1° 


Ns=80  training  vectors 
3.2  km  training  window 


note  the  narrow  clutter 
notch  of  the  3D-STAR  filter 


-40 - 1 - 1 - 

-tt.5  -0,4  -13  -0.2  -01  0  Q.1  07 

Normalized  Doppler 


Optirram 
3D  STAR 

PiE-PcppI? 


D.3  0.4  as 


L=2,  J=2  for  3D-STAR  filter  -  L=7  for  STAR  filter 


Peter  Parker 


Department  of  Electrical  &  Computer  Engineering 


Brigham  Young  University 


Provo,  UT 


84602 


BYU  Kglm 


Direct  path  JDOA=-20° 


f  Ns=80  training  vectors 
I  3.2  km  training  window 


note  the  narrow  clutter 
notch  of  the  STAR  filters 


Optimum 

STAR 

3D-STAR 

Prc-Dcppis 


-*0 - - 

-0.5  -0.4  -0.3 


-D.2  -01  0  fl.l  02 

Normalized  Doppler 


63  3.4  0.5 


L=2,  J=1  for  3D-STAR  filter  -  L=7  for  STAR  filter 
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Conclusions 

•  STAR  based  filtering  ideal  for  STAP  problems 
that  require  small  secondary  sample  support 

•  Easily  extended  to  handle  hot  or  range-varying 
clutter  models 

•  Simulations  with  realistic  circular  array  data  show 
promising  performance 

•  The  structured  nature  of  the  filters  leads  to 
computationally  efficient  algorithms 
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